Highly Eccentric Kozai Mechanism and GW Observation for Neutron Star Binaries 
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The Kozai mechanism for a hierarchical triple system could reduce the merger time of inner 
eccentric binary emitting gravitational waves (GWs), and has been qualitatively explained with the 
secular theory that is derived by averaging short-term orbital rotations. However, with the secular 
theory, the minimum value of the inner pericenter distance could be artificially blocked by the 
back-reaction of GW emission. Compared with traditional predictions, the actual evolution of an 
eccentric inner binary could be accompanied by (i) a higher characteristic frequency of the pulse-like 
GWs around its pericenter passages, and (ii) a larger residual eccentricity at its final inspiral phase. 
These findings would be important for GW astronomy with the forthcoming advanced detectors. 



I. INTRODUCTION 

Today, large-scale laser interferometers are under de- 
velopment to attain a world-wide network of second- 
generation GW detectors [1|. Their overall sensitivities 
will be improved by a factor of ~ 10, with drastic noise 
reduction at the lower frequency regime down to ~ 10Hz 
(l|. Accordingly, understandings of basic properties of 
potential astrophysical sources have become significant, 
more than ever. 

One of the most promising targets of these detectors 
is inspiral of a neutron star binary (NSB), and in this 
paper we focus our attention to GW observation for 
NSBs. From identified samples in our Galaxy, NSBs 
are expected to have very small residual eccentricities 
(0(1CT 5 )) around 10Hz [|[|. 

Meanwhile, it has been pointed out that the Kozai 
mechanism might play important roles for compact bi- 
nary mergers [4-6] . This mechanism works for hierar- 
chical triple systems, and oscillates pericenter distances 
of inner binaries, due to exchange of angular momenta 
between the inner/outer orbits [8|. This characteristic 
feature can be qualitatively understood with the secular 
theory for which, following a perturbative method in an- 
alytical mechanics, we effectively average out short-term 
fluctuations associated with both the inner/outer orbital 
rotations [3, [l0| • Since energy loss due to GW emission 
depends strongly on pericenter distance, the Kozai mech- 
anism can largely reduce the merger time of an inner NSB 
of a triple system. This interesting possibility has been 
actively discussed mostly with the secular theory includ- 
ing the averaging operations 0-01 (see also pj)- 

In this paper, we show that, with the secular theory, 
the minimum value of the inner pericenter distance could 
be artificially blocked by the back-reaction of GW emis- 
sion, in a causally unreasonable manner. To handle evo- 
lution of a highly eccentric inner binary, we need to prop- 
erly resolve the two orbital rotations without taking their 
averages. Using both simple analytical and numerical 
methods, we show that the actual evolution of an inner 
NSB could be largely different from that predicted by 
the secular theory. For the inner NSB, this could result 



in (i) a higher characteristic frequency of the pulse-like 
GWs around the pericenter passages, (ii) a higher resid- 
ual eccentricity at the final inspiral phase, and (iii) a 
shorter merger time. All of these changes could be more 
than one order of magnitude. Our findings (i) and (ii) 
are significant for the advance detectors and their data 
analyses. While quantitative evaluation for the merger 
rate requires detailed astronomical assumptions and is 
beyond scope of this paper, the last one (iii) indicates a 
higher merger rate for NSBs of triples in star clusters Q . 
This is because the outer third body would be frequently 
perturbed there. 

In this paper, we only discuss relativistic effects for 
hierarchical triples, but tidal effects around planets also 
depend strongly on orbital distance |12| . For extra-solar 
planetary systems {e.g. Hot Jupiters [13, LLJ]), an inves- 
tigation similar to this work would be worth considering. 



II. SECULAR THEORY 

We study evolution of a hierarchical triple system of 
point masses too, toi and m.2- We basically use the ge- 
ometrical unit with G = c = (mo + toi + itpz) = 1. The 
inner binary is composed by mo and toi, and we de- 
note its semimajor axis by oi. In the next section, we 
also introduce astrophysical units, considering to-o-toi as 
a NSB. For the outer third body m,2, we define its semi- 
major axis a2, relative to the mass center of the inner 
binary (total mass Mi = Too + toi). Likewise, we use 
the labels j = 1 and 2 for the inner and outer orbital 
elements (e.g. e\ for the inner eccentricity), and assume 
hierarchical orbital configurations with a = a\ja<i <C 1. 

First, we briefly discuss long-term secular evolution of 
the triple system in Newtonian dynamics, following the 
approach developed by von Zeipel 0. By suitably us- 
ing canonical transformations, we effectively average the 
short-term fluctuations associated with both the inner 
and outer mean anomalies l\ and I2 (the instantaneous 
angular positions of the inner/outer point masses [12j). 
The relevant Hamiltonian after the averaging operations, 
can be evaluated perturbatively with the expansion pa- 
rameter a <C 1. The leading order (quadrupole) term 



H 'qd = 
H qd 



0(a 2 ) is given by 



= C qd [(2 + 3e 2 )(l - 36> 2 ) - 15e?(l - 6> 2 ) COS2W!] 
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with C q d = m^rri\m 2 oi 1 j\Y%M\a 2 {\ 
gument of the inner pericenter wi [12|. Here we define 
9 = cos I with the relative inclination angle / between the 
inner/outer orbits. We denote the next order (octupole) 
term by H' oc {= 0(a 3 )). For our secular analysis of the 
inner binary, we keep up to this term for the gravita- 
tional perturbation externally induced by m 2 . But there 
exists a relation H' oc ex (too — mi), resulting in H' oc = 
for m = mi [10l | . Later we use this property to examine 
possible effects of the sub-leading terms. 

Next, we mention general relativistic corrections to 
the system, using the post-Newtonian (PN) expansion. 
The lowest order (1PN) term H[ for our hierarchical 
configuration is obtained after averaging the inner mean 
anomaly l± as [J-[a, @1 
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Here we introduced an auxiliary parameter Xipn f° r later 
discussions (xipn — 1 unless otherwise stated). At this 
stage, our effective Hamiltonian H' c for the secular evo- 
lution is given by 



H 'c - H 'qd 



H„, 



H' 

lpn ; 



(3) 



and the system is conservative (thus putting the sub- 
script "c" above) |3-Q. Using canonical equations and 
transformations of variables, we have e.g. 
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= 30C, 
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and (dai/dt) c = 



(l-r)sin2wi+O.T., (5) 



Here we defined G\ = 



tootoi [ai(l — ef)/Mi] and put O.T. for terms orig- 
inating from H' oc 0, ||| . 

The triple system becomes dissipative at the 2.5PN 
order, due to emission of GWs. Given our hierarchical 
configuration, the dissipation predominantly works for 
the inner binary, and we include its effects only for a\ 
and ei , using standard formulae for isolated eccentric bi- 
naries [ll|. Combining these with the conservative con- 
tributions, we can write down the final expressions for 
the secular evolution such as dwi/dt = (duii/dt) c , 
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with strong dependencies on 1 — ei, and da-ijdt = [J, 
6]. These secular equations have been widely used for 
analyzing long-term evolutions of relativistic hierarchical 
triple systems. 



III. NUMERICAL RESULTS 

In this section, we numerically discuss the Kozai mech- 
anism for relativistic hierarchical triples, first using the 
secular equations and then directly integrating the PN 
equations for three-body systems. While a triple system 
has many parameters, we fix most of them to concisely 
explain our new findings. 

In our geometrical unit, we fix the masses at M\ = 0.2, 
TO2 = 0.8, and the initial orbital parameters at a\ — 
3.57 x 10 5 , a 2 = 60ai = 2.14 x 10 7 = a 2l {i.e. initially 
a = 1/60), ei = 0.2 and e 2 — 0.6. We also set the 
initial angular variables at u>i — 7r/2 and tti = u} 2 = 
(Oi: the longitude of the inner ascending node [12]). For 
the secular theory, the remaining parameter is the initial 
inclination /;. We explore the regime Ii ~ 90° for which 
an inner binary can merge in a short time (also preferable 
for costly direct calculations). 

For actual astrophysical system, we presume that the 
inner binary is a NSB with their total mass M\ — 
2.8Mq. Then the initial axes correspond to a\ = 0.05 
AU and a 2 — a 2: i = 3AU. Below, instead of the di- 
rect time variable t, we use the effective outer rotation 
cycles A^2 = t/P 2 i defined with the initial orbital pe- 
riod P 2 i = 27ra 2i (corresponding to 1.38yr). The pri- 
mary GW frequency of a quasi-circular inner binary be- 
comes 10Hz (Mower end of the advanced detectors) at 
»1 = »lcr = 34.6. 

Since observed NSBs have nearly equal masses (with 
relative difference of < 7% [2J), we mainly set toi = m 2 = 
0.1 in the geometrical unit. For an isolated binary with 
a semimajor axis a = 0.05AU and masses too = toi = 
IAMq, the merger time due to GW emission becomes 
1.0 x 10 10 yr even for e = 0.7. 

As mentioned earlier, the octupole term H' oc van- 
ishes for too = mi. In order to safely estimate its po- 
tential effects, we subsidiary examine an excessive case 
(m , mi) = (0.11,0.09). 



A. Results with the Secular Theory 

As an example for predictions of the secular theory, 
in Fig.l, we provide the inner orbital elements a± and 
r v \ = <2i(l — ei) (pericenter distance), as functions of the 
outer cycles N 2 . Their ratio r v \ja\ is identical to (1 — e\). 
The basic parameters for this calculation are given in the 
caption. 

The inner binary merges at N 2m — 1209 that is consid- 
erably smaller than the cycles N 2m = O(10 10 ~ n ) for iso- 
lated binaries with moderate initial eccentricities [5j, |6| . 
Due to the Kozai mechanism, the inner eccentricity e\ 
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FIG. 1: Evolution of the inner semimajor axis ai and peri- 
center distance r p i = (1 — ei)ai. These results are obtained 
with the traditional secular theory. We set mi = m,2 = 0.1, 
ms = 0.8 with initial inclination U — 91°, and initial eccen- 
tricities ei = 0.2 and ei = 0.6. The inner binary merges at 
the outer cycles N2 m — 1209. 



oscillates with the minimum distances of r v \ ~ 300. The 
energy of the inner binary is radiated mostly around these 
epochs. As discussed in the literature [g, Q, the oscil- 
lation amplitude of e\ decreases gradually due to the 
1PN apsidal precession (the last term in Eq.Q), and, 
at N 2 > 1000, the inner elements evolve, as if an isolated 
binary. The binary becomes nearly circular at the final 
phase close to the merger. The two elements ei(<C 1) 
and a,\ there are related to their earlier values a\ e and 

4 (~ 1) as ei ~ (425/304) 145/242 (oi/2oi e (l - e le )) 19/12 
At the critical separation ai = a\ cr , the residual 
eccentricity becomes e\ cr = 5.3 x 10~ 3 . 

In Fig. 2, using the symbols on the solid lines, we 
show the duration A^m and the residual eccentricity e\ cr 
at a\ = ai cr for Ij ~ 90°. The results (circles) for 
(mo, mi) = (0.1,0.1) are similar to those (triangles) for 
(to , mi) = (0.11,0.09). Therefore, for the present pa- 
rameters, the octupole term plays a minor role, and the 
perturbative expansion itself is effective for the secular 
theory (see also [T3. fl5JV 



B. Direct Three-Body Calculations 

Now we move to direct three-body calculations. We 
use the PN equations of motions for spinless three-body 
systems [lGJj , and handle the three particles equivalently. 
In addition to the Newtonian and 1PN terms, we included 
the dissipative 2.5PN terms, but excluded the time con- 
suming 2PN terms (briefly discussed later). 

For numerical integration, we apply a fourth-order 
Runge-Kutta scheme with an adaptive time-step control 
[17| . We terminate our runs, when the inner semimajor 
axis decreases to a\ = a\ cr or when the separation be- 
tween too /toi becomes less than 10Mi. The later condi- 
tion reflects our perturbative (PN) treatment of nonlinear 
gravity, but no run encountered this condition. For nu- 
merical evaluation of Qj and ej (j = 1,2), we use the con- 
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FIG. 2: (a) The circles (mo = mi = 0.1) and triangles 
(mo = 0.11, mi = 0.09) represent the outer rotation cy- 
cles N2m before the mergers of the inner binaries (slightly 
displaced horizontally to prevent overlaps of symbols). The 
symbols with lines are obtained from the traditional secular 
theory. These without lines are from direct three-body cal- 
culations. For each inclination h, totally fifty runs with ran- 
dom initial mean anomalies are analyzed, and we show the 
median values (filled symbols), first (25%) and third (75%) 
quantiles (open symbols), (b); The residual eccentricity e\ cr 
of the inner binaries at the semimajor axes a± = a\ cr = 34.6 
(corresponding to the primary GW frequency of 10Hz for a 
NSB). 



secutive maximum ((l + ej)a,j) and minimum ((1 — Sj)dj) 
of the two orbits. 

For the direct calculations, we need to specify the ini- 
tial mean anomalies lj. Since three-body problem de- 
pends strongly on initial conditions, we randomly dis- 
tribute the initial mean anomalies to examine statistical 
trends of evolutions. For each initial inclination 7, and 
mass combination in Fig. 2, we made 50 runs, and evalu- 
ated their median values and first/third quantiles of the 
durations Nim and the residual eccentricities e\ cr . For 
m = toi =0.1 and 7j = 90°, we additionally made 50 
runs, including the 2PN terms, and obtained the median 
values N m — 78.3 and e\ cr = 0.136 that are close to the 
corresponding ones in Fig. 2. Therefore, for our analyses, 
the 2PN effect would not be important. 

We found that, in the direct calculations, the outer 
parameters a 2 and e 2 stay nearly at their initial values, 
in agreement with the secular theory. However, Fig. 2 
shows that the duration N 2m and residual e\ cr are totally 
different. 

To closely look at these discrepancies, we show an il- 
lustrative sample among the 50 runs for Ii = 91° and 
too = toi = 0.1. This run ended at N 2m = 52.5 with 
the residual e\ cr = 0.313 (close to the upper quantile in 
Fig. 2). If we simply use the outer cycle N 2 (as in Fig.l), 
the semimajor axis a\ comes to appear merely as a step 
function, and we cannot resolve its rapid final evolution. 
Therefore, for Fig. 3, we use the logarithm of the remain- 
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FIG. 3: Similar to Fig.l, but given from a direct three-body 
calculation. The inner binary takes ei cr = 0.313 at a\ = aicr, 
and merges at the outer cycle Afem = 52.5. For the horizontal 
axis, we use the remaining outer cycle AN = N2 m —N2 before 
the inner merger. 



ing cycles AN = A^m — A^ before the merger. 

We can see that, up to AN — 0(0. 1), the axis a\ is 
nearly a constant, but the pericenter distance r v \ has 
modulation of period ~ 1. This reflects the outer eccen- 
tric motion characterized by I2, and does not exist in the 
secular theory for which we average out the outer body, 
roughly speaking, as a continuous ring. 

At the turning point AA~ ~ 4x 1CP 2 , the quantity 1— e\ 
takes a minimum value, corresponding to r p \ = 38, and 
the axis a\ starts to decrease rapidly from ax = 3.0 x 10 5 . 
The distance r v \ stays nearly at a constant for a while. 
The inner binary evolves almost independently on the 
outer body 7712, and emits pulse-like GWs around the 
pericenter passages. This waveform has a characteris- 
tic frequency (Mi/rh) 1 / 2 /n ~ 10Hz that is ~ 30 times 
higher than the counterpart in Fig.l. This change would 
be encouraging for ground-based GW observation. 



IV. DISCUSSIONS 

Actually, the large reduction of r p \ sheds light on 
an important issue of the traditional secular theory for 
highly eccentric orbits. Here we discuss this from the 
viewpoint of GW emissions. For an orbit with 1 — e± <C 
1, GW emission is dominated at the pericenter pas- 
sages, and the radiated energy there is given as SE ~ 
— 857r(momi) 2 M 1 / /(12^r [ ), depending strongly on 
r p i [18| • Therefore, in order to lose a fraction Y of the 
total energy — momi/2ai, the close approach r v \ should 
be 
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with our model mass parameters, a fiducial value Y — 
0.2, and a\ at the turning point. 

Meanwhile, for the secular theory, we can estimate the 
local minimum of 1 — e\ from Eq.© r|. Neglecting the 
octupole terms, we obtain the minimum pericenter dis- 
tance r pl:min = ai(l - ei) m in as 
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(9) 



with the factor X = sin 2 I\ sin2wi| < 1. 

Therefore, even with the conservative setting X = 1, 
the minimum distance r p x jm i n of the secular theory is 
much larger than Eq. ([5]) , and only a small fraction Y <C 1 
of energy can be lost during the single passage. 

Now we carefully look at Eq.0. The fluctuations of 
the inner angular momentum (ex yja\(l — e 2 )) due to 
the outer body is mostly induced when the inner sep- 
aration is 0{a\) (considering the torque from the tidal 
force) and the radiation reaction does not work virtu- 
ally. In other words, aside from probability, no physical 
mechanism prohibits a nearly radial orbit (1 — ei ^C 1) 
in the separation regime. However, Eq.([7]) indicates that 
the minimum value of 1 — e\ is blocked by the radiation 
reaction, and the underlying analysis falls into causally 
unreasonable structure. 

To properly handle the evolution, we need to suffi- 
ciently resolve the inner orbital motion without making 
its average. Here the concrete position of the outer body 
would be also important to estimate the fluctuations in- 
duced to the small inner angular momentum |19| . If or- 
bital elements (e.g. 1 — e±) evolve more slowly than the 
two orbital periods, the traditional double averages might 
be efficient, but we need careful studies for highly eccen- 
tric binaries that could be critical to GW astronomy. 

Fig. 2 shows that the residual eccentricity e\ cr could be 
much larger than the estimation by the traditional sec- 
ular theory. This is closely related to the difference of 
the pericenter distances. For a quasi-circular binary, the 
residual eccentricity could be probed through the associ- 
ated phase modulation of GWs Q • For a NSB detectable 
with advanced detectors at SNR~ 15, the resolution of 
the residual value e\ cr (at 10Hz) would be Ae\ cr ~ 0.01 
Q . Interestingly, this is just between the two predictions 
in Fig. 2. While this figure is given for specific sets of 
parameters, our results would be intriguing for the up- 
coming GW detectors. 
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